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For  a  number  of  years  the  pulse  sounding  method  of  Breit  and  Tuve'-  -^ 
has  been  one  of  the  most  important  tools  for  investigating  the  structure  of 
the  ionosphere.      In  this  method  a  pulse  modulated  electromagnetic  wave  vath 
an  appropriate  carrier  frequency  is  transmitted  vertically  from  the  ground, 
toward  the  upper  atmosphere,   and  the  time  delay  of  the  reflected  pulse  as  a 
function  of  the  carrier  frequency  is  measured.      The  analysis  of  data  obtained 
by  the  method  of  Breit  and  Tuve  has  provided  much  of  the  conceptual  picture 
of  the  ionosphere  which  investigators  have  formed.     However,   from  time  to  time 
some  doubt  has  been  expressed  in  the  literature  about  the  validity  of 
certain  assumptions  used  to  relate  pulse  sounding  data  to  the  ionosphere's 
structure.     In  fact,   it  appears  likely  that  in  many  cases  the  mathematical 
analysis  ordinarily  applied  is  too  crude  even  to  justify  certain  conclusions 
which  are  almost  qualitative  in  nature  such  as  the  number  of  distinct  layers 
present   or  the  heights  of  maximum  electron  densities  during  a  given  experiment. 
Some  indication  of  the  pitfalls  ^ich  may  occur  in  the  usual  analysis,   as  well 
as  the  possibility  of  improving  the  results  will  be  discussed  here. 

This  discussion  will  be  limited  to  what  is   probably  the  simplest  model 
of  electrcOTiagnetic  wave  propagation  in  the  ionosphere  for  which  the  free 
electron  density  is   assumed  to  vary  with  height.     It  will  be  assumed  that  the 
variation  is  strictly  one-dimensional  so  that   the  ionized  medium  is  plane 
stratified.     In  addition,   the  earth's  magnetic  field  and  losses  in  the  propagating 
wave  due  to  collisions  of  the  free  electrons  will  be  neglected.     Under  these 
assumptions  the  effective  index  of  refraction  of  the  medium  has  the  f orm '-^-' *  l-^-' 
(1)  n(x,k)  -   (1-  V(x)/k^)^/2  ^ 
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where  V(x)  is  a  function  cf  the  height  x  and  is  proportional  to  the  free 
electron  density  while  k  is  the  wave  number  in  free  space: 
(2)  k  -  co/c 

in  terns  of  the  circular  frequency  co  and  the  free  space  velocity  of  light  c. 
The  field  component  u(x,k)  of  a  c.w.   propagating  in  the  ionized  medium  satisfies 
the  differential  equation 


(3) 


d^u/dx^   +  . 


k^-  V(x) 


.  u  ■  0  . 


A  standard  method  for  analysing  the  time  delay  data  obtained  from  a 
pulse  sounding  is  based  on  the  fact  that  the  time  delay  7'(k)  can  be  expressed 
roughly  by  means  of  the  relation 


(U) 


r(k) 


-  2  f  "  dx/Vg  , 


where  v    is  the  group  velocity  of  the  pulse  and  x.    is  the  so-called  reflection 
height,   at  which  the  index  of  refraction  vanishes.     Usually  this  expression  is 
given  in  terms  of  the  virtual  height  h'(k),   defined  to  be  half  the  time  delay 
multiplied  by  the  free  space  velocity  of  light.      In  terms  of  the  virtual  height 
there  is  a  relation  equivalent  to    (U): 


(5)  h'(k)     -     k  |k'-  V(x) 

o 


-1/2 


dx. 


where  x.  is  the  value  of  x  at  which  V(x)  is  equal  to  k  .  After  h'(k)  is 
given  by  the  experimental  data, the  integral  equation  (5)  can  be  solved  bv  a 
standard  method  for  the  function  x(V),  called  the  true  height t^-'.  The 
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desired  function  V(x)  describing  the  electron  density  distribution  is 
effectively  determined,  of  courae,  by  the  true  height  x(V). 

The  time  delay  T(lc)  is  only  approximately  defined  from  the  experimental 
point  of  view  because  in  a  sounding  when  the  reflected  pulse  returns  it  is 
distorted  bv  the  dispersive  effects  inherent  in  the  index  of  refraction  n(x,k). 
It  is  possible,  however,  to  define  the  time  delay  precisely^''-' ,  For  this 
purpose  consider  the  complex  reflection  coefficient  r(k)  at  the  ground 
corresponding  to  waves  propagating  according  to  the  differential  equation. (3). 
If  the  phase  of  r(k)  is  9(k)  so  that 

r(k)  •  |r(k)|exp  .ie(k)l  , 
then  the  virtual  height  is  given  by 
(6)  h'(k)  -  d©/dk. 

The  time  delay  and  thus  the  virtual  height  can  also  be  defined  equivalently 
in  terms  of  the  pulse,  even  when  the  dispersion  is  very  pronounced'-  -^ .  The 
pulse  time  delay  ^"(k)  is  the  time  delay  of  the  centroid  with  respect  to  time 
of  the  reflected  instantaneous  power  relative  to  the  centroid  of  the 
instantaneous  power  of  the  incident  oulse  used  in  the  sounding, provided 

the  incident  pulse  is  very  long. 

Another  source  of  difficulty  occurs  in  the  method  associated  with  the 
integral  equation  (5)  in  that  the  standard  solution  can  be  obtained  only  when 
V(x)  is  a  monotonically  increasing  function  of  x.  Thus,  the  nethod  applies 
only  to  a  height  in  the  ionosphere  up  to  which  the  electron  density  is  always 
increasing.  In  particular,  no  accurate  information  can  be  deduced  about  the 
various  layers  which  may  occur  nor  even  about  the  maximum  electron  density  of 
a  single  layer,  for  a  breakdown  in  the  approximation  will  generally  take  p3a  ce 
somewhat  before  the  height  of  maximum  density  is  reached. 
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In  principle,   the  reconstruction  of  V(x)  can  be  carried  out  exactly 
from  a  complete  knowledge  of  the  reflection  coefficient  r(k)  at  all  frequencies'^  ■' 
One  can  hardly  expect  this  much  data  to  be  available,   however.     Hence  there 
arises  the  question  of  ;just  how  much  information  about  V(x)  can  be  squeezed  out 
of  the  data  vrtiich  is  available.     Unfortunately,   at  the  present  time  it  seems 
extremely  difficult  to  find  a  ccHnplete  answer  to  this  question.     Some  idea 
can  be  obtained  about  what  cannot  be  asserted,  though,  by  considering  special 
cases.     It  can  be  seen,  for  example,   that  there  is  not  necessarily  a  one-to-one 
correspondence  between  the  local  maxima  of  the  virtual  height  versus  frequency 
curve  and  the  layers  occurring  in  the  ionized  medioim.     In  fact,   examples  exist 
with  a  single  maximum  in  the  virtual  hei^t  versus  frequency  curve  corresponding 
to  an  electron  density  distribution  function  V(x)  having  an  infinite  number 
of  local  maxima. 

At  frequencies  well  below  the  ultimate  critical  frequency  of  the 

ionosphere  the  virtual  height   provides  essentially  all  the  information  it  is 

possible  to  obtain  about  the  reflection  coefficient,   whose  modulus  will  then 

be  practically  one.     Consider  now  the  case  where  there  are  tw3    layers  at 

the  bottom  of  the  ionized  medium  such  that  the  electron  density  of  the  higher 

layer  is  much  greater  than  that  of  the  lower.     Let   |f-[  be  the  modulus 

and  6      the  phase  of  the  reflection  coefficient  of  the  lower  layer  by  itself, 
''l 

Let    |r_|   be  the  modulus  and  ©      the  phase  of  the  reflection  coefficient  of  the 
2  r_ 

upper  layer  by  itself.     Let  ©.     be  the  phase  of  the  transmission  coefficient 

■^1 
of  the  lower  layer  by  itself.     The  reflection  coefficient  of  the  combination 

of  the  two  layers  is 


(7)       r  -  ||r^|exp(i  ©^  )  +  \r^\   exp[i(e^  *  2  9^    ]| 

iHr.llr^l  exp[l(2e,^  .  e^^  -  e^^)] 


-1 
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If   |r   I   -  1,    (7)  takes  the  form 

(8)  r  -  iexpdiZf)  .  Jexpri(2©^     +  6^  )\ 

where  (/(  is  the  phase  of  the  complex  number  of  unit  modulus! 

(9)  exp(i(2f)  -  ||rj   exp[i(0^^-  29^^-  e^^)].lU|rJexp[i(2*,^*  0^^-  9^^  ^ij 

If  the  frequency  is  well  above  the  critical  frequency  of  the  lower 
layer  but  below  that  of  the  upper  layer,  |r-|  will  be  very  small  while  IrpI 
will  still  be  nearly  one.  In  this  case 


(10 )       r  '^  exp 


i(2©.  +  e  )1, 
-   *1   ^2  J 


One  may  interpret  the  derivative  «ith  respect  to  frequency  of  the  phase 
of  a  transmission  coefficient  according  to  the  physical  idea  of  the  pulse  group 
velocity.     The  phase  derivative  is  the  increment  of  transit  time  contributed 
by  the  medium  associated  with  the  given  transmission  coefficient.     It  follows 
from   (10)  that  in  order  to  find  the  virtual  height  appropriate  to  the  upper 
layer  by  itself     after  a  measurement  is  made  of  the  time  delay  of  a  pulse 
reflected  from  the  combinea  layers     one  must  correct  this  measurement  by 
subtracting  twice  the  incranent  of  transit  time  of  the  pulse  through  the  lower 
layer.     This  simple  correction  only  works,   of  course,   when  the  frequency  is 
far  enough  above  the  critical  frequency  of  the  lower  layer  so  that  the  quantity 
|r^ I   is  nearly  zero.     Otherwise  a  further  correction  must  be  made  involving 
the  quantity  0  ±n  (8)  and   (9). 

A.I.  Manning'--'   has  pointed  out  the  need  for  correcting  the  virtual 
height  raeaeuremoHtj  Moreover  he  has  shown  that  the  result  is  not  unique  when 
the  standard  analysis  based  on  the  integral  equation  (5)  is  used  to  obtain  the 
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information  about  the  lower  layer  needed  for  the  correction.  This  is  due  to 
the  fact  that  one  cannot  determine  from  the  usual  analysis  the  nature  of  the 
lower  layer  beyond  the  height  of  its  maximum  electron  density. 
From  the  fact  that 


(11)      e  ^  lim   J   p-  V(x)] 


)ll/2d,-  ' 


"dx-kx 

X->00 


>  t 


where  9  is  the  phase  of  the  transmission  coefficient  of  an  ionized  medium 
which  begins  at  x  =  0,  it  can  be  seen  in  general  that  the  corresponding  transit 
time     increment  is  approximately 


(12)  Ay 


/CO  0 

V(x)dx/2k  c, 


for  frequencies  well  above  the  greatest  critical  frequency  of  the  medium. 
It  follows  from  (11)  and  (l2)  that  the  correction  indicated  in  expression  (10) 
can  be  obtained  from  a  knowledge  of  the  total  number  of  free  electrons  in  the 
lower  layer. 

Clearly,  the  problem  of  obtaining  the  total  number  of  free  electrons  in 
the  lower  layer  must  be  approached  by  methods  other  than  the  standard  one  based 
on  the  integral  equation  (5).  One  possibility  that  may  lead  to  the  desired 
result  depends  on  the  approximation  of  a  reflection  coefficient  by  a  rational 
function  of  the  wave  number  k.  If  such  an  approximation  can  be  made,  a  simple 
expression  for  the  transmission  coefficient  can  be  obtained.  For  example,  if 
the  reflection  coefficient  r(k)  is  a  rational  function  of  k  with  simple  poles 
X  ,  V  «  1, ...,n  in  the  lower  half  of  the  complex  k  plane,  and  if  ^   ,   v-1, .,,,n 
are  the  roots  in  the  upper  half  of  the  complex  k  plane  of  the  equation 
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(13)    1  -  r(k)r(-k)  -  0, 

then  the  corresponding  transmission  coefficient  t(k)  is  given  by 

(lU)    t(ic)  "   J  J(k+A:  )/(k-x  ) 
v»l 

The  best  approximation  of  r(k)  by  a  rational  function  of  k  evidently 
must  be  an  approximation  based  on  the  virtual  height,  i.e.,  the  phase  of  r(k), 
in  a  range  of  frequencies  which  does  not  have  to  extend  too  far  above  the 
critical  frequency  of  the  lower  layer. 

It  should  be  mentioned  that  a  closed  form  expression  for  V(x)  can  be 

obtained  when  r(k)  is  a  rational  function  of  k,  although  the  expression  may 

[9] 
be  quite  complicated  in  some  cases.  More  generally,  the  problem  of  obtaining 

V(x)  from  an  arbitrary  reflection  coefficient  can  be  reduced  to  solving  a 

Fredholm  integral  equation*-  ^  ,     To  carry  out  the  determination  of  V(x),  first 

define 

00 


(15)    R(x+t)  -  5^  I   jr(k)exp[-ik(x-^t)]^  dk 


00 

Next  solve  the  integral   equation 

X 

(16)  R(t+y)K(x,y)dy  +  K(x,t)   ♦  R(xH)  -  0,     t  <  x 

-'-00 

for  the  function  K(x,t)  for  every  x  and  all  values  of  t  in  the  range 
-00  <  t  <  X.     The  function  V(x)  is  then  given  by  the  equation 

(17)  V(x)   -  2  dK(x,x)/dx. 
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The  Fourier  transform  of  (15)  in  many  cases  can  be  evaluated  by  the 
method  of  stationary  phase.  If  this  is  done,  the  kernel  R(x+t)  of  (16)  will 
depend  on  the  virtual  height  ami  its  derivatives  with  respect  to  the  wavs 
number  k.  In  principle,  then,  one  can  obtain  the  function  V(x)  frcrni  the 
virtual  height,  whether  or  not  the  electron  density  is  a  monotonic  function 
of  height. 

If  a  method  such  as  this  one  is  to  be  used  oroperly  to  analyze  the 
free  electron  distribution  of  an  ionized  gas,  it  is  probably  necessary  to  improve 
the  measurement  technique  for  obtaining  virtual  height  data.  To  some  extent 
this  can  be  done  by  using  a  test  pulse  which  is  as  long  as  possible  rather 
than  a  short  one  which  is  now  customarily  used  in  pulse  sounding  experiments'--'. 
A  further  check  on  the  data  might  be  provided  by  measuring  the  amount  of  the 
dispersion  of  the  reflected  pulse  at  each  frequency,  particularly  in  the 
neighborhoods  of  the  critical  frequencies  where  the  spread  due  to  dispersion 
is  sometimes  quite  large.  Since  the  dispersive  distortion  of  the  reflected 
pulse  is  an  effect  which  depends  entirely  on  the  reflection  coefficient,  a 
measurement  of  this  distortion  should  provide  additional  information  about 
the  reflection  coefficient. 
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